ApJ accepted 

Preprint typeset using L^T^ style emulateapj v. 04/03/99 



1 WGA J1226.9+3332 : A HIGH REDSHIFT CLUSTER DISCOVERED BY CHANDRA 

I. CAGNONI'•^ M. ELVIS^ D.-W. KIM^ p. MAZZ0TTA^-\ J.-S. HUANG^ 

AND 

A. Celotti' 

' SISSA, Via Beirut 4 -34138, Trieste, Italy 
^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

3 ESA fellow 



O 
O 

(N 

in 



> 

o 
o 
o 

:^ 

Oh 

6 



X 



ApJ accepted 



ABSTRACT 



We report the detection of IWGA J1226.9+3332 as an arcminute scale extended X-ray source with the Chandra 
X-ray Observatory. The Chandra observation and R and K band imaging strongly support the identification of 
IWGA J1226. 9+3332 as a high redshift cluster of galaxies, most probably at z = 0.85 ± 0.15, with an inferred 



temperature kT = keV and an unabsorbed luminosity (in a r = 120 aperture) of 1 .3; 



+0.16 ^ in4 
0.14 ^ 



erg S-' (0.5- 



10 keV). This indication of redshift is also supported by the K and R band imaging, and is in agreement with the 
spectroscopic redshift of 0.89 found by Ebeling et al. (2001). The surface brightness profile is consistent with a (3- 
model with /3 = 0.770 ± 0.025, r,. = ( 1 8 . 1 ± 0.9)" (corresponding to 101 ± 5 kpc at z = 0.89), and 5(0) = 1 .02 ± 0.08 
counts arcsec"^ . IWGA J1226. 9+3332 was selected as an extreme X-ray loud source with Fx/Fy > 60; this 
selection method, thanks to the large area sampled, seems to be a highly efficient method for finding luminous 
high z clusters of galaxies. 



Subject headings: galaxies: high-redshift — 
(1WGAJ1226.9+3332)- 
(1WGAJ1226.9+3332) 

1. INTRODUCTION 



galaxies: clusters: general - 
- X-rays: galaxies: clusters - 



- galaxies: clusters: individual 

- X-rays: individual 



Clusters of galaxies are good tracers of the large scale struc- 
ture of the matter distribution in the universe. The standard 
models of structure formation predict that the cluster distribu- 
tion and evolution is fully determined by the spectrum of pri- 
mordial perturbations and cosmological parameters ilo and A 
(e.g. Press & Schechter 1974; and later works), thus observa- 
tions of high redshift clusters constrain these parameters (e.g., 
Oukbir & Blanchard 1992). Moreover X-ray measurements of 
high redshift clusters of galaxies can place strong constraints 
on the thermodynamic evolution of the intracluster medium 
(ICM). For example the Luminosity-Temperature (Lx - T) re- 
lation at different redshifts probes the interrelated evolution of 
the cluster baryon mass and the total mass (e.g. Kaiser 1991; 
Evrard & Henry 1991; and later works). Furthermore X-ray 
and observations of the Sunyaev-Zeldovich effect of a sample 
of objects at different z may be used to obtain an independent 
estimate of Hq (for a review see e.g., Birkinshaw 1999) 

Given this potential wealth of information in the past few 
years a great effort has been made to search for high redshift 
clusters of galaxies (e.g. Rosati et al. 1998; Vikhlinin et 
al. 1998). Among the different methods. X-ray surveys have 
played the most important role, but no method allowed the iden- 
tification of more than a handful of clusters at z > 0.8. 
The finding of clusters from the X-ray data was complicated 
mainly by the low spatial and/or spectral resolution of the pre- 
vious X-ray missions. Einstein and ROSAT marked an impor- 
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tant step in the study of clusters but thanks to sub-arcsecond 
Chandra spatial resolution it is now possible to distinguish eas- 
ily between point sources and the more diffuse X-ray emission 
from clusters at any redshift. 

We presentthe Chandra observation of IWGA J1226. 9+3332 
This source is one of 16 peculiar ROSAT PSPC 
sources selected for their extremely high X-ray to opti- 
cal flux ratio (Cagnoni et al. 2000; Cagnoni et al., in 
prep.). IWGA J1226.9+3332 is a bright (Fo.i-2.4keV > 10"'^ 
ergcm-^s"') WGACAT (White, Giommi & Angelini 1994) 
source with blank fields, i.e. no optical counterparts on the 
Palomar Observatory Sky Survey to 0=21.5. The extreme 
Fx/Fy ratio that follows from this is incompatible with all 
major and common classes of extragalactic sources, including 
normal quasars, AGNs, normal galaxies and nearby clusters of 
galaxies (Maccacaro et al. 1988). Possibilities for the nature of 
these 'blanks' (Cagnoni et al. in prep.) include: 

(a) Quasar-2s, i.e. high luminosity, high redshift heavily ob- 
scured quasars, the bright analogs of Seyfert 2s; (b) Low Mass 
Seyfert-2s, that is AGNs powered by a low mass obscured black 
hole (i.e. obscured Narrow Line Sy 1); (c) AGNs with no big 
blue bump, e.g. ADAFs; (d) Isolated Neutron Stars undergoing 
Bondi accretion from the ISM (Madau & Blaes 1994); (e) 7- 
ray burst X-ray afterglows, (f) failed clusters, in which a large 
overdensity of matter has collapsed but has not formed galaxies 
(Tucker, Tananbaum & Remillard 1995); and, most relevant to 
this paper, (g) high redshift clusters of galaxies. 




Here we present strong evidence that IWGA J1226. 9+3332 
is indeed a high redshift cluster. 

We will use Ho=75 km s"' Mpc"' and qo = 0.5; errors in the 
paper represent la confidence levels, unless explicitly stated 
otherwise. 

2. CTAM)/M OBSERVATIONS 

IWGA J1226. 9+3332 was one of two 'blanks' for which we 
obtained Cycle 1 Chandra (Weisskopf et al., 1996) observing 
time. It was observed in ACIS-S configuration (Garmire et al., 
in prep.) with the backside illuminated S3 chip for a total useful 
exposure time of 9832.3 s. The data was cleaned as described 
in Markevitch et al. (2000a) and for the spectral and image 
analysis we used the latest available ACIS background dataset 
(Markevitch et al. 2000b). 

2.1. Spatial analysis 

To maximize the signal-to-noise ratio, we extracted the im- 
age in the 0.5-5 keV band (Fig. la and b). Along with a num- 
ber of faint point-like sources, IWGA J1226. 9+3332 is clearly 
seen as an extended source on an arcminute scale. It shows 
azimuthal symmetry, except for a possible excess to the north- 
west, ^ 40" from the cluster center. 

After subtracting the background, from the same regions in a 
normalized ACIS background map (Markevitch et al. 2000b), 
we extracted the X-ray surface brightness profile (Fig. 2) in con- 
centric annular regions centered on the X-ray emission peak 
and chosen in order to have 20 counts per annulus. The profile 
appears to be smooth without any obvious central excess related 
to a cooling flow (Figs. 1 and 2). 

We fitted the surface brightness profile with a standard (3 
model ^ (Cavaliere & Fusco-Femiano 1976) using the Sherpa 
(Siemiginowska et al., in prep.) modeling and fitting tool 
from the CXC analysis package CIAO 2.0 (Elvis et al., in 

25(r) = 5(0)[1 + (r/r,.)-]"^''+'/^ 

smaller binning, e.g. as in Fig. 3, leads to similar results 
'*Normalizing the background map using an area of IWGA J1226. 9+3332 observation without sources (7.56% lower background) consistent results are obtained: 
kT=11.94keVandz=0.87 



prep.). We obtained best fit values of /3 = 0.770 ± 0.025, 
r,. = (18.1 ±0.9)" (corresponding to 101 ±5 kpc at z = 0.89), 
and 5(0) = 1.02 ±0.08 counts arcsec"^ with a of 27.5 for 
26 degrees of freedom (d.o.f.). The excess to the north-west of 
the cluster is also visible in the radial profile; a drop of the sur- 
face brightness at ^ 40" is present in the radial profile for the 
north-west sector. Similar features in the surface brightness ra- 
dial profiles were detected by Chandra in nearby clusters (e.g. 
Markevitch etal. 1999;2000a;2001, Vikhlinin et al. 2001;Maz- 
zotta et al. 2001) and interpreted as signs of subclump motion. 

2.2. Spectral analysis 

We extracted an overall spectrum in a circle with r = 120" 
in the 0.5-10 keV band in PI channels, corrected for the gain 
difference between the different regions of the CCD. The spec- 
trum (Fig. 3) contains ^ 1100 net counts and we binned it in 
order to have 100 counts per bin^. Both the effective area file 
(ARF) and the redistribution matrix (RMF) where computed by 
weighting each position dependent ARF's and RMF's by the X- 
ray brightness. We fitted the spectrum in the 0.5-10 keV range 
with an absorbed Raymond-Smith model (Raymond & Smith, 
1977) using Sherpa. The source redshift is treated as unknown 
and, because of the the low statistics, no iron line or other line 
complex features are expected to be visible. 

In order to get an estimate of the temperature, we fixed 
the equivalent hydrogen column density to the Galactic value 
(Nh = 1 .38 X 10^° cm-^ Stark et al. 1992), the metal abundance 
to 0.3 times the Solar value, and we draw confidence levels for 
the T and z (Fig. 4). We find a best fit value of kT = 10.24 keV 
and z = 0.85 (x^ is 18.68 for 20 d.o.f.)'*. However, as shown in 
Fig. 4, these values are not well constrained. 

While submitting this paper we found that the same object 
had been independently identified as a cluster of galaxies in 
the WARPS survey (Ebeling et al. 2001) and observed for 
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Sunyaev-Zel'dovich effect measurement by Joy et al. (2001). 
The cluster spectroscopic redshift measured by Ebeling et al. 
(2001) is z = 0.89, which is within the errors of our esti- 
mate based on Chandra and optical/IR constraints (see Sec- 
tion 4). Fixing the redshift to this value we obtain a temper- 
ature of kT = 10.47!^ keV (kT=12.07 keV using a normal- 
ized background), which is in good agreement with that ob- 
tained from the Sunyaev-Zel'dovich measurement by Joy et al. 
(2001) (kT = 10.0!?;^ keV). The absorbed 0.5-2.0 keV flux^ in a 
r = 120" aperture is (3.0 ± 0.3) x lO"'-' erg cm"^ s~' , consistent 
with the Ebeling et al. (2001) PSPC measurement; the 0.5-10.0 
keV absorbed flux in the same aperture (Table 1) is (8.03:^^11) x 
10"'^ erg cm~^ s~^ For z=0.89 the unabsorbed bolometric and 
0.5-2.0 keV band luminosities are L = (2.2 ± 0.2) x lO'*^ erg s"^ 
and Lz(o.5-2.o) = (4.4 ± 0.5) x lO"'^ erg s~^ respectively. 

3. OPTICAL, INFRARED AND RADIO OBSERVATIONS 

We obtained an R-band image of the IWGA J1226.9+3332 
field on Feb. 2, 1997 using the SAO 1.2 m telescope on 
Mt.Hopkins. A R= 20.4 ± 0.2 galaxy is detected less than 1" 
from the X-ray centroid (Table 1). 

A K-band image of the field was obtained using NSFCam at 
the NASA IRTF on Jan. 31, 2000. The same galaxy is also 
seen, but is much brighter at K=15.5 (Fig. 5) isophotal magni- 
tude. The K-band magnitude within the core radius of the X-ray 
source is K=15.4, implying that little large scale IR emission is 
present. The R-K color of this object is 5.1 with an estimated 
uncertainty of 0.3 mag due to the poor image quality in both K- 
and R-band. 

Comparing R and K-band images, it is clear that there are 
many more objects detected in the K- than in the R-band. We 
detect 23 objects with K< 19.5 mag in the 1 .5' x 1 .5' image us- 
ing the source extractor software (SExtractor, Bertin & Amouts 
1996) versus only 4 at R. The star/galaxy classification was 
carried out morphologically using the Kron radius (Bertin & 
Arnouts 1996), and 3 objects are classified as pointhke. These 
objects are also detected in the R-band and have bluer colors 
(R-K<4.2) than the extended objects; we argue they are stars. 
We should not have to worry about the star contamination of the 
sample since there are very few stars with K>16 (e.g. Glaze- 
brook et al. 1995). The rest of the morphologically extended 
objects are galaxies. These, including the bright one located 
at the X-ray centroid, have very similar red colors (4.5 <R- 
K< 6.5), implying that they likely have similar redshifts and be- 



long to a cluster. A high-z cluster, CIG J0848H-4453, at z=1.27 
was detected in a near-lR field survey using similar color crite- 
ria in an over-dense region (Stanford et al 1997). The apparent 
asymmetry of the galaxy distribution may be an artifact of the 
longer exposure time at the center of the K-band mosaic image 
which is offset from the X-ray centroid. 

The galaxy surface density in the K-band image is clearly 
high (Fig. 5). We generate K-band galaxy number counts in our 
field in the range 15 < A" < 19.5, and compare them with those 
obtained from the near-lR field survey in this range (Gardner et 
al. 1993, Saracco et al. 1997, Huang et al. 1997, Minezaki et 
al. 1998, Huang et al. 2000) with their coverage ranging from 
200 arcmin'^ to 10 deg^. The number counts obtained from our 
image are substantially higher (> lOcr, a factor of 100 at K=19) 
than those from the field surveys. Such a large excess cannot 
be due to Poisson statistics or magnitude errors. 

Both the high density and the similar (extremely red) colors 
for these galaxies thus imply that they are likely to be members 
of a cluster. 

The FIRST survey detected a faint 3.61 ± 0. 1 8 mJy source at 
1.4 GHz (Becker, White & Helfand, 1995) close to the center 
of the X-ray emission (Table 1). The radio source is pointlike 
(FWHM< 0.91", < 16 kpc at z = 0.85) and has a luminosity 
Lr(z = 0.85) ~ 6.6 X 10^4 W Hz'^, compatible with low lumi- 
nosity radio-loud AGNs (e.g. Zirbel & Baum, 1995). 

4. DISCUSSION AND CONCLUSION 

Chandra has shown that IWGA J1226.9+3332 is an ex- 
tended X-ray source, with a hard, high temperature ther- 
mal spectrum. Optical and IR imaging has shown that 
IWGA J1226. 9+3332 has faint optical/lR counterparts. A 
cluster of galaxies is the only known type of object that could 
fit such a description. Moreover the K-band image shows a 
strong excess of galaxies compared with field counts at K> 17 
(Fig. 5) around IWGA J1226.9-I-3332 . Physical clustering is 
the only possible explanation. Several lines of argument go on 
to suggest that it is a high redshift cluster of galaxies. Below 
we list these arguments and try to constrain its temperature and 
redshift. These results are summarized in Figure 4. 
(1) The Chandra X-ray profile is well fitted by a /? model with 
(3=0.11, a value in agreement with a typical relaxed cluster 
(e.g. Jones & Forman, 1999 and references therein). More- 
over, if the cluster redshift is 0.7 <z< 1-2, then the observed 
angular core radius, r, = (18.1 ±0.9)" (101 ± 5 kpc at z = 0.89 
with the assumed cosmology), corresponds to a linear size of 



'The main component in the flux error is the uncertainty on the cluster temperature (quoted in the text); using the ztlu values on T we derive an error of +11% 
-12% on the flux. Other contributions come from the normalization chosen for the background (±9%), from the fraction of counts lost using the r = 120" aperture 
assuming a /3-profile (±5%) and from the Poisson error on the counts (±3%) 
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Fig. 3. — Chandra ACIS Spectrum and residuals to a Raymond-Smith model with kT = 10.23 keV and z = 0.89. The spectrum was binned, for 
display purposes, to obtain a minimum of 30 counts per bin. 

Fig. 4. — The contours of the fit with a Raymond-Smith plasma at 1, 2 and 3 a. The three blue solid lines represent the luminosity-temperature 
relation (Markevitch 1998) for qo = 0.5 and the three blue dashed lines for qo = 0. The red lines represent the limits on the redshift obtain from the 
R-K color of the brightest object assuming it is a first ranked elliptical (Coleman, Wu & Weedman 1980), while the dotted lines are the limits from 
the magnitude of the first ranked elliptical (see text). The big dot corresponds to the best fit to the Chand ra spectrum. 



90 < r < 150 kpc for any value of (130 < r < 220 kpc for any 
value of Vl for //o=50 km s~' Mpc"'), values consistent with a 
typical relaxed cluster. 

(2) It is well known that clusters of galaxies follow a well- 
defined Luminosity-Temperature relation (e.g. Markevitch 
1998 and reference therein). Recently it has been shown that 
the local Lx-T relation does not evolve (or is consistent with 
little evolution) with redshift up to z « 0.8 (see e.g. Wu et al. 
1999; Delia Ceca et al 2000; Fairley et al. 2000 and reference 
therein). Fig. 4 shows that the Lx-T relation (L = AT^ where 
T6 = T/6 keV, a = 2.02 ± 0.40 and A = (1 .71 ±0.21) x lO^"* 
erg s"') from Markevitch (1998) requires a 3(t lower limit of 
r > 4 keV and z>0.4, while the Icr limits require T> 7.5 keV 
and z > 0.65 if no evolution is assumed. The best fit value of 
0.85 obtained from Chandra spectrum (dot in Figure 4) is con- 
sistent with the unevolving Lx-T relation. 

(3) K-band imaging has shown that IWGA J1226.9+3332 has a 
galaxy at K= 15.5. If this source is a first ranked cluster elliptical 
with Mk = -26.7 ± 0.5 (Collins & Mann 1998), then, assuming 
negligible K-correction, qo = 0.5 and , Ho=50 km s"' Mpc"' as 
in Collins & Mann (1998), it has z = 0.68!j] (dotted lines in 
Fig. 4). 

(4) the color of the K-band galaxies is extremely red, with none 
bluer than R-K^ 4.5 and a maximum R-K= 6.5. Only a un- 
evolving elliptical galaxy at 0.7 <z< 1.5 or a Sbc galaxy at 
z> 1.1 can have such a red color (Coleman, Woo & Weed- 
man, 1980). Using the more accurate R-K^ 5.1 ±0.3 of 
the first ranked elliptical we can restrict the redshift range to 
0.75 < z < 1 .0 (red lines in Fig. 4). 

Using all these constraints we conclude that 
IWGA J1226. 9+3332 is a distant cluster of galaxies with a 
most probable redshift of 0.85 ±0.15 and not smaller than 
z = 0.65. 

This gives kT = 10^3 keV and implies an X-ray luminosity, 
determined in a r = 120" aperture, of Lx(o.5-iokeV) = l-3^o|4 x 



10"*^ erg s"', corresponding, for z = 0.85, to a bolometric L = 
(2 ± 0.2) X 10"^^ erg s"'. Our estimated redshift is similar to the 
spectroscopic z = 0.89 found by Ebeling et al. (2001). 

The blank field X-ray source IWGA J1226.9+3332 is thus 
a highly luminous and massive high redshift cluster and a use- 
ful source to determine the evolution of the clusters X-ray lu- 
minosity function (e.g. Rosati et al., 1998). Since models in 
the direction of low il universe (with or without cosmological 
constant) (e.g. Henry 2000, Borgani & Guzzo 2001) predict 
a higher density of high redshift clusters compared to high fl 
models, finding hot high redshift clusters has a strong leverage 
on cosmological models. 

Since such high luminosity, high redshift clusters should be 
rare, the relative ease with which this discovery was made is 
potentially of great significance. The search for high Fx/Fy 
sources ('blanks'), sampling a large area of the sky, is an ef- 
ficient method of finding very luminous high redshift clusters; 
serendipitous typical flux limited surveys can find (and found) 
plenty of z > 0.6 clusters, but they are inefficient at finding 
such luminous clusters because they cover relatively small ar- 
eas (^ 100 deg^). This methodology is a useful complement to 
serendipitous flux limited surveys. 
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TABLE 1 

IWGA J1226. 9+3332 OBSERVATIONS 



Energy Band 


Instrument 


Date 


Exposure 


Coordinates 


Offset" 


Count Rate 


Flux (units) 








(ks) 


(J2000) 


(") 


(Counts s"') 




X-ray 


Chandra 


Jul 31 2000 


9832.3 


12 26 58.2-1-33 32 48.28 


0.0 


0.107 ±0.006* 


8.03 X 10-" (ergcm-^s-'O 


R-band 


SAO 48in 


Feb 02 1997 


900 


12 26 58.2 -1-33 32 48.7 


0.87 




20.4 ±0.2 (Mag.) 


K-band 


IRTF 


Jan 31 2000 


1200 


12 26 58 -H33 32 48^ 


2.0 




15.5 (Mag.) 


Radio 


FIRST 






12 26 58.19-H33 32 48.61 


0.79 




3.61 ±0.18 (mjyat 1.4 GHz) 



" Offset from Chandra position 

* [0.5-10 keV] in a circle with r=120" computed from the spectral model within sherpa '' Due to the lack of bright stars in the K-band image, we obtained an estimate 
of the position using the objects in common with the R image. 



